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Mirror dark matter offers a framework to explain the existing dark matter direct 
detection experiments. Here we confront this theory with the most recent exper- 
imental data, paying attention to the various known systematic uncertainties, in 
quenching factor, detector resolution, galactic rotational velocity and velocity dis- 
persion. We perform a detailed analysis of the DAMA and CoGeNT experiments 
assuming a negligible channeling fraction and find that the data can be fully ex- 
plained within the mirror dark matter framework. We also show that the mirror 
^ ■ dark matter candidate can explain recent data from the CDMS/Ge, EdelweissII and 

CRESSTII experiments and we point out ways in which the theory can be further 
tested in the near future. 
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1 Introduction 



The field of dark matter direct detection has blossomed in recent times, with ex- 
citing positive signals from DAMAfl, 2, 3], CoGeNT[4], as well as interesting hints 
from CDMS/Ge[5] and CDMS electron scattering [6]. Very recently, more exciting 
evidence for the direct detection of dark matter has arisen from the CRESSTII[7] 
and EdelweissII[8] experiments. 

Mirror dark matter has emerged as a simple predictive framework which can 
explain all of the direct detection experiments [9, 10, 11, 12, 13]. The purpose of this 
article is to provide a comprehensive update of the experimental status of the mirror 
dark matter candidate, paying particular attention to the various known systematic 
uncertainties, in quenching factor, detector resolution, galactic rotational velocity 
and velocity dispersion. 

Recall, mirror dark matter posits that the inferred dark matter in the Uni- 
verse arises from a hidden sector which is an exact copy of the standard model 
sector [14] (for a review see ref. [15]). That is, a spectrum of dark matter particles of 
known masses are predicted: e', H', He', O', Fe', ... (with m e > = m e , mu> = itih, etc). 
The galactic halo is then presumed to be composed predominately of a spherically 
distributed self interacting mirror particle plasma comprising these particles [16]. 
In addition to gravity, ordinary and mirror particles interact with each other via 
(renormalizable) photon-mirror photon kinetic mixing[14, 17]: 

fornix — > (-Q 

where F^ (F' ) is the ordinary (mirror) U(l) gauge boson field strength tensor. 
This interaction enables mirror charged particles to couple to ordinary photons 
with electric charge q = ee and thus allows mirror particles to elastically scatter 
off ordinary particles. This means that mirror dark matter can be probed in dark 
matter direct detection experiments. It turns out that this simple predictive theory 
can explain the DAMA annual modulation signal, the CoGeNT low energy excess 
as well as hints from CDMS, Edelweiss and CRESSTII consistently with the null 
results of other experiments. 

The outline of this paper is as follows. In section 2 we provide a brief review 
of the mirror dark matter theory. In section 3 we provide some necessary technical 
details: cross section and halo distribution which are characteristic of mirror dark 
matter. In section 4 (5), we examine the implications of the most recent DAMA 
(CoGeNT) data for the mirror dark matter theory. These experiments are sensitive 
to dark matter particles heavier than around 10 GeV which makes them excellent 
probes of the dominant mirror metal component of the galactic halo, A'. We show 
that these experiments can be simultaneously explained and lead to a measurement 
of the parameters: ey 7 ^ and tua 1 both of which are consistent with the theoretical 
expectations of e ~ 10~ 9 (from galactic halo energy balance) and A' ~ O' =>- itla' ~ 
16m p (from analogy with the ordinary matter sector). In section 5 we also show 
that the DAMA and CoGeNT signals are consistent with the results of the other 



1 



experiments including the null results of XENON100 and CDMS/Si. In section 6 we 
examine the constraints on e' scattering from the DAMA absolute rate. We show 
that these constraints when combined with the DAMA and CoGeNT data suggest a 
halo mirror metal proportion £4/ ~ 1CT 2 . In section 7 we examine recent data from 
the CDMS/Ge, EdelweissII and CRESSTII experiments. CDMS/Ge and Edelweiss 
are excellent probes of the anticipated Fe' component, and the data are consistent 
with a Fe' component with mass fraction: ^Fe'l^A' ~ 10~ 2 . We also point out that 
the CRESSTII experiment is potentially sensitive to both A' and Fe' components 
and their recently announced low energy excess can be explained by A' and Fe' 
interactions. In section 8 we draw our conclusions. 

2 A brief review of mirror dark matter 

Mirror dark matter conjectures that the inferred dark matter in the Universe arises 
from a hidden sector which is an exact copy of the standard model sector. That is, 
the standard model of particle physics is extended: 



Such a theory can be theoretically well motivated from symmetry considerations if 
left and right handed chiral fields are interchanged in the extra sector. In this way 
space-time parity symmetry and in fact the full Poincare group can be realized as 
an unbroken symmetry of nature, and for this reason we refer to the particles in the 
extra sector as mirror particles. The standard model extended with a mirror sector 
was first studied in ref . [14] and shown to be a phenomenologically consistent renor- 
malizable extension of the standard model. The concept, though, has a long history 
dating back prior to the advent of the standard model of particle inter actions [18]. 
For a review and more complete list of references see ref. [15]. 

If we include all interaction terms consistent with renormalizability and the sym- 
metries of the theory then we must add to the Lagrangian a U(l) kinetic mixing 
interaction[17, 14] and Higgs - mirror Higgs quartic coupling[14]: 



where F^ u (F' ) is the ordinary (mirror) U(l) gauge boson field strength tensor 
and (<f)') is the electroweak Higgs (mirror Higgs) field. The most general Higgs 
potential, including the quartic Higgs mixing term (above) was studied in ref. [14] 
and shown to have the vacuum (0) = (0') for a large range of parameters. With 
this vacuum, the mirror symmetry is unbroken and consequently the masses of the 
mirror particles are all identical to their ordinary matter counterparts. 

In this framework, dark matter is comprised of a spectrum of stable massive 
mirror particles: e',H',He',0'... etc, with masses m e > = m e ,mH> = %, "Ihv 
in if, etc. To explain the rotation curves in spiral galaxies, the dark matter needs 
to be roughly spherically distributed in galactic halos. Given the upper limit on 



£ = CsAi(e, u, d, 7, ...) + C SM (e', u', d', 7', ...) . 



(2) 



C mx = \F^F^ + A0t#V , 



(3) 
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compact star sized objects (MACHOs) in the halo from microlensing observations, 
roughly fmacho ~ 0.2 — 0.3 depending on the assumptions [19], we then expect the 
mirror particles to be distributed predominately as a hot gaseous spherical halo 
surrounding the collapsed disk of ordinary matter [16] 2 . 

Observations of colliding clusters, such as the bullet cluster [20] indicate that dark 
matter does not have self interactions on galaxy cluster scales. This suggests that 
the gaseous mirror dark matter component is confined to galactic halos (c.f. [21]). 
Gravity and the mirror particle self interactions may well be sufficient to achieve 
this. 

A dissipative dark matter candidate like mirror matter can only survive in an 
extended spherical distribution in galaxies without collapsing if there is a substantial 
heating mechanism to replace the energy lost due to radiative cooling. In fact, 
ordinary supernova can plausibly supply the required heating if the photon and 
mirror photon are kinetically mixed with e ~ 10~ 9 [16] 3 . For kinetic mixing of this 
magnitude about half of the total energy emitted in ordinary Type II Supernova 
explosions (~ 3 x 10 53 erg) will be in the form of light mirror particles „ T , e>± i l') 
originating from kinetic mixing induced plasmon decay into e' + e'~ in the supernova 
core [23]. Given the observed rate of Supernova's in our galaxy of about 1 per century, 
this implies a heating of the halo (principally due to the e /=t component), of around: 

Lllt-in ~^ 3 x 10 53 erg 1 ~ 10 44 erg/s, for the Milky Way . (4) 

It turns out that this matches (to within uncertainties) the energy lost from the halo 
due to radiative cooling[16]: 

L h e ner 9 y-out = A / n\ Anr 2 dr ~ 10 44 erg/s, for the Milky Way. (5) 

In other words, a gaseous mirror particle halo can potentially survive without col- 
lapsing because the energy lost due to dissipative interactions can be replaced by the 
energy from ordinary supernova explosions. Presumably there are feedback mecha- 
nisms which maintain this balance. For example if L^ l ° rgy _ out > L^ at _ in then the 
halo would contract which in turns increases the gravitational pull on the ordinary 
matter component. This compression of the ordinary matter component should in- 
crease ordinary star formation rates, thereby increasing Lf^ at _ in until the energy is 
balanced. In this way the ordinary supernova rate might be dynamically adjusted so 
that the halo is stabilized. Extending these ideas to galaxies beyond the Milky Way, 
the hypothesized connection between Supernova rates and dark matter distribution 



2 Naturally a MACHO subcomponent consisting of mirror white dwarfs, mirror neutron stars 
etc are also expected and can be probed by microlensing observations. Since most of the stellar 
mass is ejected as gas in the explosions producing these stellar remnants, it is plausible that the 
MACHO mass fraction can satisfy the observational limit of fmacho ~ 0.2 — 0.3. 

3 A mirror sector with such kinetic mixing is consistent with all known laboratory astrophysical 
and cosmological constraints [22]. 
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might ultimately provide a dynamical justification for the empirical Tully-Fischer 
and Faber- Jackson relations. 

The spherically distributed mirror particle plasma is likely to be far too hot 
for much mirror star formation to occur at the present epoch. Thus, we do not 
expect significant heating of the ordinary matter sector from mirror Supernovas at 
the present time. However during the first billion years or so the situation might 
have been the reverse. That is, the early stages of galaxies may have witnessed 
rapid mirror star formation and evolution and little ordinary star formation, due 
ultimately to the effects of asymmetric initial conditions in the Early Universe. In 
particular, for T' <C T in the early Universe (required to achieve successful big bang 
nucleosynthesis and large scale structure formation) 4 and e ~ 10~ 9 the primordial 
mirror He' abundance is expected to be relatively high, Yue ~ 0.9 [24]. With such 
initial conditions the evolution rate of stars is dramatically increased by several 
orders of magnitude [25]. In other words, we surmise that the required asymmetric 
evolution of the ordinary and mirror matter components in galaxies, originates in 
the complex interactions between the ordinary and mirror components (such as 
the energy transfer to the mirror sector from ordinary supernova's) together with 
asymmetric initial conditions in the early Universe. 

Since mirror charged particles have an electric charge induced via the kinetic mix- 
ing, mirror nuclei can interact with ordinary nuclei via spin independent Rutherford 
elastic scattering. It turns out that the positive results of DAMA and CoGeNT 
can be explained by such scattering from a putative ~ O' component if e ~ 1(U 9 . 
This provides important experimental evidence in favour of the mirror dark matter 
candidate, which we now examine in detail. 

3 Interaction cross-section and galactic mirror dark 
matter distribution 

The interaction rate in an experiment depends on the cross-section, da/dEpt, and 
halo velocity distribution, f(v). The photon-mirror photon kinetic mixing enables a 
mirror nucleus [with mass and atomic numbers A', Z' and velocity v] to elastically 
scatter with an ordinary nucleus [presumed at rest with mass and atomic numbers 
A, Z] . In fact the cross-section is just of the standard Rutherford form corresponding 
to a particle of electric charge Ze scattering off a particle of electric charge eZ'e. The 
cross-section can be expressed in terms of the recoil energy of the ordinary nucleus, 
E R [9]: 

da A 
dE~R = £W 1 

4 See ref.[26] for further discussions about early Universe cosmology with mirror dark matter. 



4 



(6) 



where 



A = Fi(gr A )Fl,(gr>) , (7) 

and Fx(qrx) {X = A, A') are the form factors which take into account the finite size 
of the nuclei and mirror nuclei. [The quantity q = (2i71aEr) 1 ^ 2 is the momentum 
transfer and rx is the effective nuclear radius] 5 . A simple analytic expression for the 
form factor, which we adopt in our numerical work, is the one given by Helm[27, 28]: 

Fx(qrx) = 3^1-^ , ( 8 ) 
qr x 

with rx = 1.14X 1 / 3 fm, s = 0.9 fm and ji is the spherical Bessel function of index 
1. 

The halo mirror particles are presumed to form a self interacting plasma at an 
isothermal temperature T. This means that the halo distribution function is given 
by a Maxwellian distribution: 

fi(v) = e-^ 2 /T 

= e -" a /«gW ? ( 9 ) 

where the index % labels the particle type [i = e' , H' , He' ,0' , Fe' , ...] and v 2 [i] = 
2T/mi. The dynamics of such a mirror particle plasma has been investigated 
previously[16, 9], where it was found that the condition of hydrostatic equilibrium 
implied that the temperature of the plasma satisfied: 

T * \mv 2 rot , (10) 

where rh = J2 n i m i/J2 n i is the mean mass of the particles in the plasma, and 
v rot rs 254 ± 16 km/s is the galactic rotational velocity of the Milky Way [29]. The 
velocity dispersion of the particles in the mirror particle halo evidently depends on 
the particular particle species and satisfies: 

« = ■ (11) 

rrii 

Note that if mi ^> m, then v 2 [i] <^ v 2 , ot . Consequently heavy mirror nuclei have 
their velocities (and hence energies) relative to the Earth boosted by the Earth's 
(mean) rotational velocity around the galactic center, v rot . This allows a mirror 
nuclei in the 'oxygen' mass range ~ m ~ 15 GeV to provide a significant annual 
modulation signal in the energy region probed by DAMA (Ef> > 2 keVee) which, as 
we will see, has the right properties to fully account for the data presented by the 
DAMA collaboration [1, 2, 3]. 



5 Unless otherwise specified, we use natural units, U = c = 1 throughout. 
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According to the above considerations, in order to calculate the velocity dis- 
persion, Vq[i], we need to estimate the mean mass of the mirror particles in the 
plasma. The plasma is expected to be completely ionized since it turns out that 
the temperature of the plasma is T \ keV. We start by making the simplifying 
assumption that the mirror metal component of the plasma is dominated by a single 
element, A'. Under this assumption, the plasma consists of e', H', He' and A'. It is 
straightforward to estimate m: 

m 1 

where & = ; — — : is the halo mass fraction of species i, m v is the 

n H im p +n He imHe+n A im A i ^ ' P 

proton mass and A' is the mass number. Thus, combining Eq.(ll) and Eq.(12) we 
have 

2 / V rot 
Uo ^] = TM9 57 7 t {I 3yT • ( 13 ) 
A [2- ^He'+U'KAf ~ 2)\ 

If we vary £a' between and 1, then we obtain a lower and upper limit for t>o[^4']: 

i < m. < _j_ ( i4) 

A {2- I^He') V rot 1 + — 

Mirror BBN studies [24] indicate that the primordial Helium mass fraction, Yhv, is 
relatively high, with Yud — 0.9, which is quite unlike the case of ordinary matter. 
However, like the ordinary matter sector, the primordial value for £4' is expected [24] 
to be small £4/ <^ 1. That is, heavy mirror elements are anticipated to be synthesised 
in mirror stars. If the net A' production from mirror stars remains subdominant, 
i.e. £a> <^ 1, then we expect that vq[A!\ to be given by the lower limit in Eq.(14) 
with ^He' = Y He i — 0.9. The situation where £4' ~ 1 corresponds to extremely 
efficient mirror star formation and evolution, and is a priori possible. We will there- 
fore consider both limiting cases for t>o[^4']- It turns out, though, that our results 
are relatively insensitive to the possible variation of i>o[^4'] given by Eq.(14) simply 
because v [A']/v rot « 1. 



4 The DAMA experiment 

The DAMA experiments [1, 2, 3] employ large mass scintillation sodium iodide de- 
tectors operating in the Gran Sasso National Laboratory. These experiments were 
initially operating with a target mass of around 100 kg and since 2003 with a target 
of ~ 250 kg. These experiments have been running for more than 13 years and now 
have a cumulative exposure of 1.17 ton x year. Importantly, the DAMA experiments 
have consistently observed a positive dark matter signal, with statistical significance 
of around 8.9 a C.L.[3]. 
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The DAMA experiments utilize the annual modulation signature, which provides 
a "smoking gun" signal for dark matter. The idea[30] is very simple. The interaction 
rate must vary periodically since it depends on the Earth's velocity, v E , which 
modulates due to the Earth's motion around the Sun. That is, 



where v & = v rot + 12 km/s is the sun's velocity with respect to the galactic halo and 
Ave — 15 km/s, oo = 2n/T (T = 1 year) with to = 152.5 days (from astronomical 
data). Importantly the phase and period of the modulation are both predicted! 6 
This gives a strong systematic check on their results. Such an annual modulation 
has been found in the 2-6 keVee recoil energy region at the 8.9<r confidence level, 
with T, t measured to be [3]: 



Clearly, both the period and phase are consistent with the theoretical expectations 
of halo dark matter. There are no known systematic effects which could produce the 
modulation of the signal seen and thus it is reasonable to believe that the DAMA 
experiments have detected dark matter. 

Mirror dark matter explains the DAMA annual modulation signal via kinetic 
mixing induced elastic (Rutherford) scattering of the dominant mirror metal com- 
ponent, A', off target nuclei. [The H' and He' components are too light to give a 
signal above the DAMA energy threshold] . We leave it up to the experimental data 
to determine m^/, although our best 'theoretical' guess would be A' ~ O' given that 
O is the dominant metal in the ordinary matter sector. The differential interaction 
rate is given by: 



where N T is the number of target atoms per kg of detector (we must sum over 
Na and I interactions separately), k = (tcVq[A'}) 3 ^ 2 is the Maxwellian distribution 
normalization factor and riA> = Pdm£,A' / m A' is the number density of the mirror 
nuclei A' at the Earth's location (we take pd m = 0.3 GeV/cm 3 ). Here v is the 
velocity of the halo particles relative to the Earth and is the velocity of the 
Earth relative to the galactic halo. Note that the lower velocity limit, v m i n (Ef>), is 

6 Deviations from to = 152.5 days are possible if there is bulk halo rotation. However, a large 
deviation from the expected value of to = 152.5 days would be difficult to reconcile with the inferred 
approximate spherical distribution of the dark matter in the galactic halo. 




(15) 



T = 0.999 ± 0.002 year 
to = 146 ± 7 day. 



(16) 




(17) 
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given by the kinematic relation: 



{m A + m A ,) 2 E R ^ 



\ 2m A m\, 



The velocity integral in Eq.(17), 



/ = r i^i^v (19) 

J\v\>v min (E R ) fe|V| 

can easily be evaluated in terms of error functions assuming a Maxwellian distribu- 
tion: f A ,(v,v E )/k = (irv%[A'})- 3 / 2 ex P (-(v + v E ) 2 /v%[A'}). In fact, 

1 = 2^] [erf{x + V) ~ 6rf{x ~ V)] ' (20) 

where 

_ v min (E R ) v E 

vo[A>] ' y -v [A'} ■ {21) 

The differential interaction rate, Eq.(17), can then be expanded in a Taylor series 
yielding a time independent part (which we subsequently denote as the 'absolute' 
rate) and time dependent modulated component: 

dR dRP dR^ , . . „ 

cosw(t-to) , (22) 



dE R dE R dE R 
with 

dR° N T n A ,\I(E R ,y ) 

dE R El 

dR 1 N T n A ,\Ay (dp 

dE~R El [djj 



(23) 



y=yo 



Here y = v Q /v [A'], Ay = Av E /v [A'] and 

e -(*-w) a + e -(*+W)) a l _ (24) 



dpi I(E R ,y ) | 1 

dy)y= yo Vo V^y v [A'} 



To compare with the measured rates we must take into account the quenching 
factor and detector resolution. We include detector resolution effects by convolving 
the rates with a Gaussian: 

dR^ 1 1 f dR° 



dE% V2^a res J dE R 



f d^ e -(E R -E^/2a^ dER j (25) 

J dER 
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where is the measured energy. The resolution is given by[31] 



v. 



res 



a 



+ P 



(26) 



E R ^E R (keVee) 



where a = 0.448 ± 0.035, (3 = (9.1 ± 5.1) x 10~ 3 . The unit of energy is the electron 
equivalent energy, keVee. For nuclear recoils, in the absence of any channeling, 
keVee = keV/gA, where qA is the quenching factor. For DAMA, q Na 0.3 while 
qi ~ 0.09. Channeled events where target atoms travel down crystal axis and planes 
have qA — 1. 

The issue of channeling has recently been re-examined in ref.[32]. It was found 
that the channeling fraction is likely to be very small (< 1%) in the energy range of 
interest in contrast with the earlier study[33] performed by the DAMA collaboration. 
It is argued that the DAMA analysis did not take into account that the scattered 
atoms originate from lattice sites and hence cannot be easily channeled. In light of 
these developments, we expect that the channeling fraction is indeed small - probably 
negligible. Throughout this paper, therefore, we shall generally assume that no 
channeling occurs, with the exception of figures 2 and 4 where we compute the 
allowed region of parameter space and include the channeling region as a comparison. 

The measured annual modulation amplitudes for the 1.17 ton-year cumulative 
exposure, S^io^, are binned into AE = 0.5 keVee energy bins and can be obtained 
from figure 6 of ref. [3]. This is to be compared with the computed annual modulation 
amplitude for mirror dark matter obtained by averaging the differential rate over 
the binned energy taking into account the resolution and quenching factors: 7 



We consider the energy range 2-8 keVee, separated into 12 bins of width 0.5 keVee, 
which encompasses the 2-6 keVee DAMA signal region. Varying the parameters 
ui a', ey / 6?" around the best fit, we can obtain the DAMA allowed region 8 . There are 
a number of systematic uncertainties which can be included in the analysis and we 
have examined the following: a) considering i>o[/4/] within its expected limits given 
by Eq.(14), b) varying the quenching factors for Iodine and sodium by ±20%, i.e. 
taking q^ a = 0.30 ±0.06 and qi = 0.09 ±0.02, c) varying the detector resolution over 

7 The DAMA data are efficiency corrected so there is no need to include the detection efficiency 
in Eq.(27). 

8 For the purposes of the fit, we analytically continue the mass number, A' , to non-integer values, 
with Z' = A'/2. Since the realistic case will involve a spectrum of elements, the effective mass can 
be non-integer 




(27) 



It is convenient to define a x 2 quantity: 




(28) 
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its 2a uncertainty and d) varying v ro t = 254 ±32 km/s, that is around ±2a from its 
estimated value. The variation of quenching factor and resolution were taken into 
account by minimizing x 2 {Qii QNa, Ores-, m A', ^\f^A') ov er 20% variation of qi and q^ a , 

_2 

and over the 2a uncertainty in a res . This defines X (ra^, eV^A')- The best fit for 

_2 

DAMA has X min — 8.5 for 10 degrees of freedom. Some examples near the best fit, 
assuming no channeling occurs, are shown in figure 1. 
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Figure 1: DAMA annual modulation amplitude versus measured recoil Energy for 
the parameters: mA>/m p = 20, e^/^A' = 7.4 x 10~ 10 and v rot = 254 km/s. The 
solid (dashed) line corresponds to the lower (upper) vq\A!\ limit given in Eq.(14). 
Negligible channeling has been assumed. 

Favoured regions in the mA'^y/iA' plane can be obtained by evaluating the 

contours with X =X min +9 (roughly 99% C.L.). In figure 2 we plot the allowed 
regions for DAMA assuming v rot = 222 km/s [fig 2a], v rot = 254 km/s [fig 2b] and 
v ro t = 286 km/s [fig 2c] . Also shown in the figures are the DAMA allowed region as- 
suming that channeling occurs with fractions as originally estimated by the DAMA 
collaboration [33]. It has been emphasised recently[34] that the systematic uncer- 
tainty in q Na might be as large as q^ a = 0.30 ± 0.13 given the lack of measurements 
of the quenching factor in the low energy region. If this is the case, then we find 
that the favoured regions extend out to somewhat 10 — 15%) lower values 
than that given in figure 2. 



10 



1e-08 



< 1e-09 



1e-10 



v rot = 222 km/s 




10 



20 30 
m A ./m p 



40 



50 



Figure 2a: DAMA [99% C.L.] allowed region in the mA'i^yf^A 1 plane, assuming 
negligible channeling fraction, for v rot = 222 km/s. Solid (dashed) line corresponds 
to the lower (upper) v [A'] limit given in Eq.(14). Also shown for comparison is the 
DAMA allowed regions if channeling occurs with fractions originally estimated by 
the DAMA collaboration. 
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Figure 2b: Same as figure 2a, except with v rot = 254 km/s. 
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Figure 2c: Same as figure 2a, except with v rot = 286 km/s. 



The main features of the annual modulation spectrum predicted by A' dark 
matter can be easily understood. As discussed earlier [10], at low Er, where x(Er) <C 
y, dR}/dEn is negative. As Er increases, dB}/dER changes sign and reaches a 
maximum at the value of Er where x(Er) ps y, or equivalently, the value of Er 
where v min (E R ) = v rot . At high E R (x > y), dR 1 /dE R — >■ 0. From Eq.(18) this 
means that the position of the peak, E R eak , is given by: 

ppeak _ 2m A m 2 A , 2 

h R 7 1 \2 V rot ■ V 2 ^ 

(m A + m A ,y 



For low itla' ~ 30m p the annual modulation signal arises predominantly from A' 
scattering with Na, while for m A i ~ 30m p scattering off both Na and I contributes 
significantly to produce the signal 9 . Note that since t>o[^4'] v ro t the velocity 
distribution is so narrow that the width of the peak is dominated by the detector 
resolution. This explains the relative insensitivity of the fit to the particular v [A'] 
value which is evident in figures 1,2. 

To summarize, we see that without channeling the DAMA annual modulation 
signal can be explained for a relatively wide range of m A i values, which includes 
the region around m A i/m p ~ 16, expected if A' = O' dominates the mirror metal 



9 In the case where channeling is assumed with the fractions originally estimated by the DAMA 
collaboration [33] the annual modulation signal is dominated by interactions with / only. See ref. [10] 
for further discussion of this case. 
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component. We now turn to the positive low energy excess observed by CoGeNT, 
together with the implications from the null experiments, such as XENON100 and 
CDMS/Si. 



5 The CoGeNT experiment 

The CoGeNT experiment operating in the Soudan Underground Laboratory has 
recently presented new results in their search for light dark matter interactions [4]. 
With a low energy threshold of 0.4 keVee and a Germanium target, they have 
observed a low energy excess which is not easily explainable in terms of known 
background sources. The energy region probed by CoGeNT overlaps with the en- 
ergy region in which the DAMA collaboration have observed their impressive annual 
modulation signal and thus it is natural to interpret the CoGeNT low energy ex- 
cess in terms of A' mirror dark matter interactions. In ref.[13] we showed that the 
CoGeNT excess is compatible with mirror dark matter expectations and thus pro- 
vides a model dependent check of the DAMA signal, which we now examine in more 
detail. 

To compare with the measured event rate, we include detector resolution effects 
and overall detection efficiency: 

^1 = c 1 f ^ e -(E R -E^)V2^ dE m 

dE% n Rl ^aJ dE R H K ' 

where E R is the measured energy and a 2 = cr 2 + (2.35) 2 E R i]F with a n = 69.4 eV, 
7] = 2.96 eV and F = 0.29[4, 35]. The detection efficiency, e f (E%), was given in 
figure 3 of ref.[4], which we approximate via 

^ = WW < 31 > 

The energy is in keVee units (ionization energy). For nuclear recoils in the absence 
of any channeling, keVee = keV/q, where q ~ 0.21 is the relevant quenching factor 
in the near threshold region [34]. 

We fit the CoGeNT data in the low recoil energy range assuming A' dark matter 
and that the background is an energy independent constant, together with two 
Gaussians to account for the 65 Zn (1.1 keV) and 68 Ge (1.29 keV) L-shell electron 
capture lines. Initially fixing mA>/m p = 20 and v rot = 254 km/s, as an example, 
we find a best fit of Xmm — 14.8 for 21 — 4 = 17 degrees of freedom, with ey 7 ^ = 
6.2 x 10~ 10 (independently of whether we take the upper or lower limiting values of 
vq). This fit for CoGeNT is shown in figure 3. 
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Figure 3: Fit of the CoGeNT spectrum for m A ,/m p = 20, t^/U> = 6.2 x KT 10 and 
v rot = 254 km/s. Thick solid (dashed) line is with lower (upper) v [A'] limit given 
in Eq.(14). Also shown is the dark matter contribution to the signals (thin lines). 
Negligible channeling has been assumed. 



The shape of the spectrum is nicely fit by A' dark matter due primarily to the Er 
dependent Rutherford cross section: da/dE R oc l/E R . As further data is collected 
this E R dependence will be more stringently constrained, which will pose a more 
rigorous test of the mirror dark matter theory. 

As in the DAMA case, to account for some of the possible systematic uncertain- 
ties we vary the quenching factor by ±20% [i.e. take qc e = 0.21 ± 0.04]. We define 
a x 2 function, 

x 2 (qce, e\/&' m A') = S ( j^k ~ datai J / a i ( 32 ) 

where is the differential rate averaged over the binned energy. We take the 

Ft 

experimental errors to be purely statistical, so that <7j = V 'datai. We minimize 

_2 

X 2 { ( lGei m A'i cVCa*) over the 20% variation in qc e , which defines X {ttla', eVd')- Of 
course, we also minimize x 2 with respect to the parameters of the background model 
describing the amplitudes of the 65 Zn (1.1 keV) and 68 Ge (1.29 keV) L-shell electron 
capture lines and adjusting also the constant background component. Note that no 
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background exponential is assumed (or needed) to fit the data. The CoGeNT allowed 

_2 _2 

region in the m^ey 7 ^ plane is then defined via contours X =X min +9 (roughly 
99% C.L.). 

In figure 4 we plot the allowed regions for CoGeNT together with the DAMA 
allowed region for v rot = 222 km/s [fig 4a], v TOt = 254 km/s [fig 4b] and v rot = 286 
km/s [fig 4c] . We have allowed for a 20% systematic uncertainty in both CoGeNT 
and DAMA quenching factors. In these figures we have assumed £4/ < 1 so that 
i>o[^4'] is given by the lower limit given in Eq.(14). The case where £4' — 1 features an 
almost indistinguishable allowed region (as already illustrated in figure 2 for the case 
of DAMA), so is not given. Also shown are the corresponding 95% C.L. exclusion 
limits from CDMS/Si[36], CDMS/Ge[5] and XENON100[37] experiments. 




Figure 4a: DAMA and CoGeNT [99% C.L.] allowed regions in the m A '/m p , 
plane, assuming negligible channeling fraction, for v ro t = 222 km/s. Also shown are 
a) DAMA allowed regions if channeling occurs with fractions originally estimated 
by the DAMA collaboration and b) the exclusion limits from CDMS/Si, CDMS/Ge 
and XENON100 experiments. [The region excluded is to the right of the exclusion 
limits] . 
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Figure 4b: Same as figure 4a, except with v rot = 254 km/s. 
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Figure 4c: Same as figure 4a, except with v rot = 286 km/s. 
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In computing the exclusion limits we have allowed for a 20% systematic uncer- 
tainty in energy threshold. That is, the threshold for CDMS/Si[36] was taken to be 
8.4 keV nr rather than the quoted 7.0 keV nr , the threshold for CDMS/Ge[5] was taken 
to be 12 keV nr rather than the quoted 10 keV nr and the threshold for XENON100[37] 
was taken to be 10.4 keV nr rather than the quoted 8.7 keV nr 10 . For the CDMS/Si 
experiment, we used the quoted[36] raw exposure of 33.9 kg-days for Er < 10 keV nr , 
38.8 kg-days for 10 < En(keV nr ) < 15 and 53.5 kg-days for Er > 15 keV nr , and 
assumed a detection efficiency of 20% in the energy region near threshold. For the 
CDMS/Ge experiment we assumed the total raw exposure of 1010 kg-days and a 
detection efficiency of e/ = 0.18 + 0.007 Er in the low energy region of interest [5]. 
For the XENON100 experiment we used the quoted [37] raw exposure of 447 kg-days 
together with a detection efficiency of 0.4 in the low energy region. 

From figure 4 it is evident that the CDMS/Si experiment is the most sensitive of 
the null experiments to A' dark matter, which is due to its light target element and 
relatively low threshold. The CDMS/Si experiment constrains the allowed region to 
mA'/rrip ~ 30. Observe also that the CDMS/Si experiment seems to exclude the two 
events seen in the CDMS/Ge experiment [5] as being due to the same component 
which can explain the DAMA and CoGeNT data. It is possible, though, to interpret 
the two events seen by CDMS/Ge as a hint of a heavier ~ Fe' component. This 
interpretation requires[12] ^Fe'/^A' ~ 10~ 2 and is not excluded by CDMS/Si or any 
other experiment. Such a small Fe' component does not significantly affect the fit 
of the DAMA or CoGeNT experiments. We will examine the effects of the Fe' 
component for the higher threshold experiments in more detail in section 7. 

We have also performed a global analysis of the DAMA and CoGeNT signals. 
Fixing v rot = 254 km/s and evaluating x 2 fc> r the combined DAMA+CoGeNT data, 
we have found Xmin = 24 - 4 for 28 d.o.f. at m A ,/m p = 24, e\f^X> = 6.4 x 10~ 10 . 
Excellent fits to the combined DAMA and CoGeNT data are also obtained for 
v ro t = 222 km/s and v ro t = 286 km/s. In figure 5 we plot the favoured region of 
parameter space for the combined fit of the DAMA and CoGeNT data for three 
representative values of v rot . 

10 Exclusion limits from the XENON10 experiment are only marginally better than the 
XENON100 exclusion limit (but not better than the CDMS/Si limit). For a recent discussion 
about the calibration and other uncertainties in the XENON low energy region, see ref.[38]. 
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Figure 5a: DAMA and CoGeNT [99% C.L.] global allowed region in the mA'/rrip, ey 7 ^ 
plane, assuming negligible channeling fraction, for v rot = 222 km/s. Also shown are 
the exclusion limits from CDMS/Si, CDMS/Ge and XENON100 experiments. [The 
region excluded is to the right of the exclusion limits]. 
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Figure 5b: Same as figure 5a except that v rot = 254 km/s. 
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Figure 5c: Same as figure 5a except that v rot = 286 km/s. 

The CoGeNT experiment has been running continuously since December 2009, 
and might collect enough data to search for the annual modulation signal. In figure 
6 we show results for the annual modulation amplitude predicted for the CoGeNT 
experiment. Interestingly we see that there is a change of sign for the annual modula- 
tion amplitude at low energies Er ~ 0.5 — 0.8 keVee (depending on the parameters). 
This means that at the lowest energies CoGeNT should see more events during the 
(northern) winter/fall than the (northern) summer/spring. This might provide a 
useful means of experimentally distinguishing this dark matter theory from other 
possible explanations. 
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Figure 6: CoGeNT annual modulation amplitude versus recoil energy. For some 
representative parameters: ey^A> = 7.0 x 10~ 10 , m A ijm v = 18 (solid line), ey/^/ — 
6.4 x 10" 10 , m A ,jm v = 20 (dashed line), e^/U> = 6.0 x 10~ 10 , m A ,jm v = 22 (dotted 
line). All for v rot = 254 km/s. The figure assumes 100% detection efficiency. 



6 Electron scattering and the DAMA absolute 
rate 

In addition to nuclear recoils, electron recoils from mirror electron scattering off 
bound atomic electrons in the target volume can also occur. The e' component of 
the halo is expected to be distributed with Maxwellian velocity distribution, and 
since m e i <C rh, we expect v [e'] ^> v rot . In fact, from Eqs.(ll,12) we can estimate 
that vo[e'] is in the range: 

~ m p 

Since vo[e'] ^> v rot we can, to a good approximation, neglect v rot when dealing with 
e' — e scattering. This means that electron recoils give a negligible contribution to the 
annual modulation signal, however they can contribute significantly to the absolute 
energies below 2 keVee. Interestingly such a rise at low energies was observed 
in the CDMS electron scattering data [6] and is compatible with e ~ 10 _9 [11]. Recall 
that the DAMA experiment doesn't discriminate against electron recoils, so the 
DAMA experiment is also sensitive to electron recoils through their contribution to 
the absolute rate as we now discuss. 
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2 — %£h p i 
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m r 



1 + A! 12 



(33) 
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The DAMA experiment constrains the absolute event rate to be less than around 
1 cpd/kg/keVee in the region near threshold. We shall here examine the implications 
of this constraint for mirror dark matter. As discussed above, we expect both nuclear 
recoils and electron recoils to contribute to the absolute rate from dark matter 
interactions. The nuclear recoil contribution can easily by calculated as in Eq.(25). 
The electron recoil contribution is more complicated. An accurate treatment of 
the cross section would require knowledge of the wavefunctions of all the electrons 
in the Na and / atoms. A rough estimate of the electron scattering contribution 
can be made by [11] considering only the contribution of the loosely bound (binding 
energy less than 0.1 keV) outer shell electrons in Na and /. The number of such 
loosely bound atomic electrons is 9 for Na and 17 for /. We further approximate 
these electrons as free and at rest, and compute the elastic scattering rate on these 
electrons. Thus, within this approximation the cross section has the form given in 
Eq.(6), with A e = 2Tce 2 a 2 /m e . The predicted differential interaction rate is then: 

dR n AT m f da f e ,(v) LJ 3 



dE R J dE R k 



gN T n e i I — — — \v\d v 



= gN T n e A r ^TT^ (34) 

E R J\v\>v min (E R ) k\v\ 

where is the number of target Nal pairs per kg of detector and k = (^v 2 [e'}) 3 ^ 2 
is the Maxwellian distribution normalization factor. The quantity g = 26, is the 
number of loosely bound atomic electrons per Nal pair as we discussed above. 
Also, n e i is the halo e' number density. Assuming the halo is fully ionized it is 
straightforward to show that 



1 _ ^ He> _ idl 
2 2 



Pdm (35) 



m p 



Note that the lower velocity limit in Eq.(34), v m i n (E R ), is given by the kinematic 
relation: 



' 2Er (36) 



The velocity integral in Eq. (34) can be analytically solved leading to: 

dR „ A e / 2e 



dE R = sNTn "K{^WV (37) 

where x = v m i n /vo[e']. Finally, to compare with the experimentally measured rate 
we convolve this rate, with a Gaussian [as in Eq.(25)] to take into account the finite 
detector resolution. 

There are many potential systematic uncertainties in the absolute rate, and we 
consider the following: the uncertainty in the measured detector resolution [i.e. 
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taking a 2a variation of the measured a res given in Eq.(26)], a ~ 30% uncertainty in 
the e' — e scattering cross section, and a 0.25 keVee uncertainty in energy calibration. 

In figure 7 we give an example of the absolute rate predicted for the DAMA 
experiment showing both electron and nuclear recoil contributions separately, using 
the aforementioned systematic uncertainties to minimize the rate. Since the data 
below 2 keVee is formally below the DAMA threshold we do not attempt to fit this 
data, and must await the forthcoming DAMA upgrade which is designed to lower 
the energy threshold. The rise in event rate below 2 keVee, which is illustrated in 
figure 7, is a prediction of this model which DAMA can potentially confirm when 
they lower their energy threshold. 




4 5 
Energy [keVee] 



Figure 7: DAMA absolute rate from nuclear recoils with negligible channeling (dot- 
ted line), electron recoils (dashed line) and combined (solid line) for the parameters: 
m A ,/m p = 20, e^/U> = 5.0 x lO" 10 and e = 1.0 x 10" 9 (=► £ A , = 0.25). v rot = 222 
km/s has been assumed. 
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In figure 8, we given an example of the electron scattering rate predicted for 
the CDMS Germanium electron scattering experiment [6], for the same parame- 
ters as chosen for figure 7, together with a simple linear model for the background 
(R(back ground) = 1.9 — 0.09Er). For the Germanium experiment the resolution is 
given by [6] 

a = ^(0.293)2 + (0.056) 2 E R /keV keV (38) 
and the rate as in Eq.(37) but with g — 14[11]. 11 . 



10 



8 





1 I I I I 


v rot = 222 km/s 










1 1 1 1 1 





1 2 3 4 5 6 7 8 

Energy [keVee] 

Figure 8: CDMS/Ge absolute rate from electron recoils for the same parameters as 
figure 7, together with a simple linear background model. 



The rise in event rate seen below 2 keV in both DAMA and CDMS electron 
scattering data is very interesting, but is formally below the threshold of both of 
these experiments, and therefore needs to be confirmed by future measurements. 
Conservatively, the only limit that we can obtain is by looking at the data above 2 
keV. Demanding that the total rate in the first energy bin above threshold be less 
than 1 cpd/kg/keV, suggests an upper limit for e. This upper limit depends quite 
sensitively on the value of i>o[e'] (and hence also on v rot and £4/). If £4/ <C 1, then 

11 In ref.[ll], we assumed a cutoff at Er = 0.8 keV. However this is, in fact unjustified, and here 
we have no such cutoff except for a phenomenological cutoff at Er = 0.2 keV. Such a low energy 
cutoff is necessary due to the divergence of the cross section in the Er — > limit. 
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we find: 



e ~ 2.7 x 1(T 9 , if v rot = 222 km/s 

e ~ 1.5 x 10~ 9 , if v rot = 254 km/s 

e ~ 1.0 x 1(T 9 , if v rot = 286 km/s (39) 

If £4/ ~ 1 then we find: 

e ~ 1.1 x 1(T 9 , if v rot = 222 km/s 
e ~ 0.8 x 10~ 9 , if v rot = 254 km/s 

e ~ 0.6 x 10~ 9 , if v rot = 286 km/s. (40) 

In computing these upper limits we have allowed for some of the systematic un- 
certainties by including: the uncertainty in the measured detector resolution [i.e. 
taking a 2a variation of the measured a res given in Eq.(26)], a ~ 30% uncertainty in 
the e' — e scattering cross section, and a 0.25 keV uncertainty in energy calibration. 
It should be emphasized, though, that the systematic uncertainty can potentially 
be much larger in view of the large event rate at low energies which is smeared into 
Er ~ 2 keV by the resolution. In particular the resolution has not been measured 
at Er ~ 2 keVee and the naive extrapolation might breakdown at low energies. It is 
also possible that the resolution might fall off faster than a Gaussian in the tails of 
the distribution, which would weaken the above limits on e. Finally, astrophysical 
uncertainties in modelling the halo will add further systematic uncertainties to the 
e' — e scattering rate due to its sensitive dependence on i>o[e']. Departures from 
spherical symmetry or a rotating halo etc will lead to deviations from Eq.(10), and 
hence to v [e']. [Note though that the A' scattering rate is much less sensitive to 
uncertainties in v since v [A'] <C v rot , and thus such uncertainties will have little 
affect on our DAMA/CoGeNT fit]. Thus, given these systematic uncertainties our 
limits on e should be viewed more as a guide, than strict upper limits. 

Similar limits to the above, with the same caveats regarding potentially larger 
systematic uncertainties, can be obtained from the CDMS electron scattering data. 
In combination with our estimate of ey^T ps (7 ± 3) x 10~ 10 from the DAMA and 
CoGeNT experiments, the limits given in Eq.(40) indicate £4/ ~ 10~ 2 at v rot = 222 
km/s, with stronger bounds at higher v rot values. This suggests that the mirror 
sector may have a higher metal content than the ordinary matter sector. This is 
certainly possible, and might be due to a period of rapid mirror star formation and 
evolution during the first few billion years of the Universe (which is suspected given 
the computed high primordial Y He i ^0.9 abundance [24] which would dramatically 
speed up mirror star evolution by several orders of magnitude [25]). 



7 CDMS/Ge, EdelweissII and CRESSTII 

The CDMS/Ge and EdelweissII experiments utilize a Germanium target, and both 
of these experiments have found evidence for dark matter interactions. Due to their 
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relatively high threshold of 10 keV nr for CDMS/Ge and 20 keV nr for Edelweiss, these 
experiments are not sensitive to the the dominant A' component. The light mass and 
narrow velocity dispersion of the A' component ensure that the predicted rate for 
these experiments is much less than 1 event for their net exposures of approximately 
200 kg-days and 322 kg-days respectively. However, these experiments are sensitive 
to heavier mirror dark matter components of the halo, and provide the most sensitive 
probes of the anticipated Fe' component. 

In figure 9 we have given an example of the recoil energy spectrum predicted for 
CDMS/Ge (figure 9a) and for EdelweissII (figure 9b). The numerical work assumed 
the CDMS/Ge[5] (EdelweissII [8]) resolution was given by a res = 0.2 keV (a res = 1.0 
keV) and detection efficiency, eff = 0.18 + 0.007£ R (eff ~ 1). The value for v [Fe'] 
was obtained from Eq.(13) assuming that £4/ <C 1, i.e. 



v 2 

u rot 



We'} = , ™ - , (41) 

01 J [m Fe /m p ][2-\i Hel ] V ; 



where ^He 1 ~ 0.9. 
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Figure 9b: Edelweiss spectrum for Fe 1 dark matter with ey/£pe> — 5.0 x 10 -11 . We 
have assumed v rot = 254 km/s. 



The CDMS/Ge experiment finds two events at energies 12.3 and 15.5 keV which 
are clearly compatible with the shape of the predicted recoil energy spectrum for 
CDMS/Ge. Edelweiss finds 2 events in their acceptance region with energy just 
above threshold, which is also compatible with the shape of the predicted recoil 
energy spectrum for EdelweissII. It is therefore plausible that both of these exper- 
iments have detected Fe' dark matter. Under this assumption, we can obtain an 
estimate of the parameter: e\/^Fe' for each of these experiments: 

ey/W = (3.5l£{j) x 1CT 11 from CDMS/Ge 

eyffd = (1.4lJ;jj) x 1(T 10 from EdelweissII . (42) 

Here we have only included the statistical errors at 95% C.L. The systematic uncer- 
tainties can be quite large, given the rapidly rising event rates towards lower Er. 
For example, we find that a 20% systematic uncertainty in energy scale would lead 
to a ~ 20% uncertainty in the estimate for ey/^p e ' from CDMS/Ge and a 50% in the 
estimate for ey/^p e > from EdelweissII. Also note that the systematic uncertainties 
in the form factor begin to be quite significant for Edelweiss due to the large recoil 
energy threshold. Clearly systematic uncertainties can reconcile the two estimates 
of ey/Fe' from CDMS/Ge and EdelweissII. The XENON100 experiment, with an 
anticipated net exposure of over 1000 kg-days should be able to confirm the presence 
of a Fe' signal in the near future. 
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Combining the above estimate for ey/^Fe' with our earlier fit for ey/^A> suggests 
a £,Fe'/£,A' fraction of around ~ 10~ 2 which is plausible. Also note that such a small 
C,Fe' component does not significantly affect the fit of the DAMA and CoGeNT 
experiments. 

The CRESSTII experiment using a CaWO^ target has recently reported 32 dark 
matter candidate events, with a background of around 9 events in their signal region, 
which suggests a statistically significant low energy excess of around 23 events. The 
threshold of CRESSTII is 10 keV nr , with the excess of events reported in the oxygen 
band near threshold. The CRESSTII experiment is potentially sensitive to both the 
A' component and the Fe' component. We illustrate this in figure 10, where we give 
the predicted CRESST recoil energy spectrum for an example with parameters close 
to the DAMA/CoGeNT best fit. As this figure illustrates, the A' component is only 
important in the region near threshold, while the Fe' contribution is somewhat more 
spread out. Thus, in principle these two components can be distinguished from the 
observed energy distribution of events (so long as they both contribute significantly 
to the signal). 
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Figure 10: CRESSTII spectrum in the oxygen band for A', Fe' dark matter assuming 
v rot = 254 km/s. The dashed line is the A' contribution with parameters ey / £^7 = 
6 x 10~ 10 , niA' = 22.0. The dotted line is for the Fe' contribution with parameters 
eV^Fe' — 1-0 x 10~ 10 . The solid line is the sum of the two contributions. [100% 
detection efficiency has been assumed, and a resolution of 0.3 keV]. 
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We have found that the A' contribution has a very large annual modulation in 
the energy region above threshold. In fact, 

^ « 0.5^ . (43) 
dE% dE% 1 ; 

The large annual modulation results because the interactions arise from A' particles 
in the tail of the narrow Maxwellian distribution and therefore can be greatly affected 
by small changes in the velocity of the Earth. Thus, if A' does contribute significantly 
to the signal, an examination of events in the Ethreshold < E < Ethreshold+3 keV 
bin should show a statistically significant annual modulation with just 1-2 years of 
data. The annual modulation predicted for the Fe' events is smaller, but still quite 
large: 

dR dR , . 

— — - w 0.16—— (44) 
dE% dE% 

and could also be eventually seen provided that a significant proportion of the events 
are due to Fe! interactions. 



8 conclusion 

In conclusion, we have confronted the mirror dark matter theory with the most 
recent experimental data. We examined the DAMA experiments allowing for the 
possibility of a negligible channeling fraction and showed that under that assump- 
tion the impressive DAMA signal can be fully explained. Mirror dark matter can 
simultaneously explain the CoGeNT low energy excess, and remains compatible 
with the results of the other experiments, including interesting hints of dark matter 
detection from CDMS/Ge, EdelweissII and CRESSTII. Taking into account some 
of the possible systematic uncertainties in quenching factor, detector resolution, 
galactic rotational velocity and velocity dispersion, we have mapped out the allowed 
regions of parameter space in the rriA>,ey/^A' plane [Figures 4,5]. The net result 
is that the mirror dark matter candidate can explain all of the existing direct de- 
tection experiments, with parameters ey^T = (7 ± 3) x 1CT 10 , mA>/m p = 22 ± 8, 

U> ~ io- 2 , MU> ~ io- 2 . 

Importantly this theory will soon be more stringently tested by: a) Further data 
from DAMA: in the near future the DAMA collaboration plan to upgrade their 
experiment with the aim of lowering their energy threshold. As illustrated in figure 
1, they should see a change in sign of their modulation amplitude between 1.0-2.0 
keVee. b) The CRESSTII experiment, with light target element O' and threshold 
10 keV is potentially sensitive to both A' and Fe' components. They should find a 
rapidly falling energy spectrum with most of their events between 10 and 14 keV, 
with a very large annual modulation, c) The CDMS/Si experiment is sensitive to 
the dominant mirror metal component, A'. This experiment currently provides the 
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strongest constraint from the null experiments and limits mA'/Tn p ~ 30. Further 
data from CDMS/Si should either see a signal, or produce tighter constraints on 
vtia'- d) More data from the CoGeNT experiment should be helpful. Of particular 
note, is that mirror dark matter predicts a detectable annual modulation signal for 
CoGeNT (figure 6) with the distinctive feature that it changes in sign at energies 
around Er 0.5 — 0.8 keVee. e) Very sensitive but typically higher threshold 
experiments, such as CDMS/Ge, EdelweissII and XENON100 can potentially probe 
the heavier ~ Fe' component which should exist at some level. The two events seen 
in CDMS/Ge and in EdelweissII are consistent with this component, and suggests 
£,Fe'/£,A' ~ 10~ 2 . f) In the longer term, directional experiments will be important 
due to the low velocity dispersion of mirror dark matter [fo[^4'] <S v rot \. 12 

Mirror dark matter is not expected to show up in collider experiments through 
the photon-mirror photon kinetic mixing induced interactions. However, the Higgs 
- mirror Higgs quartic interaction [Eq.(3)] leads to modifications of the properties of 
the Higgs boson[14] which can potentially be observed at the LHC and Tevatron[39]. 
Sensitive orthopositronium studies[40], might be able to directly probe e ~ 10~ 9 and 
thus provide further tests of the mirror dark matter scenario. Such an experiment 
has been proposed recently in ref.[41]. 
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